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INTRODUCTION 


Reported  here  are  four  areas  in  which  we  have  made  significant 
progress  in  the  past  year  (1  November  1984-31  October  1985) . 

1.  TISSUE  COMPOSITION  DETERMINATION  USING  ULTRASONIC  MEASUREMENTS  AND 
MIXTURE  RULES. 

This  work  is  appended  to  this  progress  report.  It  should  be  noted 
that  the  approach  discussed  is  not  restricted  to  tissues.  We  shall  con¬ 
sider  other  compositions  (e.g.,  composite  rubbers)  in  the  near  future. 


2.  MICROPARTICLE  CHARACTERIZATION  USING  ACOUSTIC  SCATTERING. 

Mr.  Roy,  Prof.  Apfel 

The  "microparticle  characterization"  project  involves  the  develop¬ 
ment  of  an  acoustic  scattering  technique  for  determining  the  mechanical 
properties  of  small  (s5  pm)  particles.  We  scatter  tone  bursts  of  30  MHz 
center  frequency  and  2  psec  duration  off  of  individual  particles  (e.g., 
biological  cells)  as  they  pass  through  the  focal  region  of  two  confocally 

i 

positioned  acoustic  transducers.  The  scattered  echos  are  detected  at 
90°  and  180°.  With  this  information  we  employ  Rayleigh  scattering  theory 
to  calculate  the  compressibility  and  density  of  the  individual  particles, 
provided  we  have  some  a  priori  knowledge  of  the  scatterer  size. 

Instrument  development  is  progressing  smoothly.  The  transducer 
positioning  system  and  the  flow  system  (for  convecting  the  particles) 
are  up  and  running.  The  transmit/receive  switch,  which  permits  us  to 
detect  backscattered  information  using  the  same  transducer  that  generates 


the  initial  tone  burst,  is  finally  operational.  We  have  also  developed 


and  debugged  the  necessary  instrumentation  and  software  for  demodulating 
the  received  signals,  sampling  the  peak  signal  amplitude,  digitizing  the 
results,  and  storing  the  data  in  the  computer. 

At  this  point,  the  data  analysis  is  carried  out  by  the  computer. 

We  seek  to  develop  a  scheme  for  "real  time"  data  analysis  employing 
analog  circuits.  This  will  provide  instantaneous  particle  information 
which  may  be  used  to  trigger  some  type  of  downstream  activity  (such  as 
sorting).  We  also  hope  to - incorporate  a  Coulter  orifice  downstream  of 
the  confocal  region.  This  device  would  provide  size  information  for 
individual  particles,  thereby  increasing  the  accuracy  of  the  scatter¬ 
ing  calculations.  We  can  calibrate  the  system  by  using  "reference 
particles"  such  as  polystryene  spheres  or  liquid  drops  and  we  plan  to 
run  experiments  with  biological  particles  such  as  white  blood  cells. 

3.  INTERFACE  CHARACTERIZATION  USING  MODULATED  RADIATION  PRESSURE  ON 
ACOUSTICALLY  LEVITATED  DROPS. 

We  had  to  make  changes  in  the  design  of  our  levitation  cell  so 
that  surfactants  could  be  cleaned  out  readily.  (They  tend  to  hide  in 
cracks  and  crevices.)  First  we  changed  the  original  lucite  cell  to  a 
rectangular  glass  cell  (  51  *  51  *  150mm,  made  by  Vitro  Dynamics,  Inc.) 
with  a  hollow  tube  transducer  (  1.5  *  1.5inch  )  glued  on  the  bottom  of 
the  cell  by  Torr  seal.  The  major  advantages  of  the  new  system  are  that 
it  is  very  easy  to  clean  even  though  there  are  surfactants  in  it,  and 
it  allows  us  to  study  the  large  amplitude  oscillations  without  worrying 
about  cavitation  or  the  dissolving  of  the  silver  paint  on  the  transducer 


The  new  system  is  also  very  efficient.  It  required  3  to  4  volts  to 
levitate  a  hexane  drop  (density  0.66  g/crP)  in  water  with  the  old 
system,  but  it  only  requires  2  volts  using  the  new  system. 

We  now  have  begun  to  study  the  large  amplitude  oxcillations  using 
the  new  cell.  The  dynamics  of  the  vibrating  system  becomes  nonlinear 
when  the  amplitude  of  the  oscillation  gets  sufficiently  large.  Dif¬ 
ferent  modes  of  oscillations  are  coupled  together  causing  the  shift 
of  the  resonance  frequency.  And  the  drop  behaves  essentially  as  a 
soft  spring  due  to  the  intrinsic  nonlinearity  of  the  governing  equa¬ 
tions  of  motion  and  the  boundary  conditions?  i.e.,  the  resonance 
frequency  decreases  as  the  amplitude  increases.  Thus  far,  the  system 
has  been  investigated  qualitatively,  because  of  difficulties  we  have 
encountered  during  the  experiments.  In  order  to  study  the  system 
quantitatively,  we  need  to  be  able  to  measure  three  major  parameters: 
the  amplitude  of  the  oscillation,  the  pressure  amplitude  of  the  sound 
field,  and  the  resonance  frequency.  The  resonance  frequency  can  be 
measured  without  any  problem,  in  the  same  way  as  for  the  small  ampli¬ 
tude  oscillation.  To  measure  the  amplitude  of  the  oscillation,  the 
drop  should  be  in  a  fixed  position,  but  that  has  not  been  the  case. 

The  drop  tends  to  move  around  as  the  amplitude  increases.  Furthermore, 
one  can  no  longer  use  the  voltage  across  the  transducer  as  an  indicator 
of  the  intensity  of  the  sound  field,  because  the  response  of  the  trans¬ 
ducer  may  not  be  linear  for  large  working  voltages.  Hence  a  calibrated 
transducer  is  needed  to  probe  the  sound  field  directly. 

If  the  amplitude  of  the  oscillation  is  increased  further,  the  drop 
will  fission.  This  can  be  achieved  very  easily  by  adding  surfactants  to 


reduce  the  surface  tension.  We  have  observed  drop  fission  for  a  trans¬ 
ducer  ve]!a;-,p  of  15  volts,  or  so.  As  expected  for  the  quadrupole  mode, 
drop  fission  results  in  two  drops  of  equal  size. 

4.  NONPRCPAGATING  HYDRODYh’AHIC  SOLITON. 

E.  Carr  Everbach,  Prof.  R.  E.  Apfel 

The  properties  of  a  nonpropagating  hydrodynamic  soliton  (as  originally 
investigated  by  Wu  and  Rudnick)  have  been  investigated  experimentally  for 
certain  resonator  geometries.  In  particular,  solitons  in  annular  resona¬ 
tors  of  constant  width  (7.5  cm)  but  different  radii  of  curvature  have 
profiles  which  deviate  progressively  from  the  hyperbolic  secant  function 
as  the  mean  radius  of  the  annulus  is  decreased.  This  distortion  reaches 
a  limit  when  the  circumference  of  the  inner  wall  becomes  on  the  order  of 
the  soliton  width,  and  the  soliton  appears  to  separate  into  two  component 
oscillations  on  either  side  of  the  inner  wall.  These  oscillations  have 
nearly  secant  profiles  and  demonstrate  most  of  the  characteristics  observed 
for  solitons  in  rectangular  resonators.  When  the  radius  of  the  inner  wall 
is  decreased  further,  the  two  component  oscillations  retain  their  relative 
separation  distance  regardless  of  the  inner  wall  radius ,  and  form  two 
secant-profile  bumps  when  no  inner  wall  is  present.  In  this  case,  the 
annular  resonator  becomes  a  circular  tank,  but  it  is  surprising  that  the 
oscillations  do  not  correspond  to  any  of  the  normal  modes  of  the  ■fank  and, 
in  fact,  occur  at  a  drive  frequency  slightly  below  that  of  the  lowest 
normal  mode. 

Recent  efforts  have  been  devoted  to  trying  to  understand  this  sub-(0.2) 


mode,  and  to  determine  its  relationship  to  the  nonpropagating  soliton. 


Theoretical  analyses  of  both  the  soliton  and  the  normal  modes  of  the 
circular  tank  have  been  attempted,  but  due  to  the  high  nonlinearity 
inherent  in  the  phenomenon,  this  has  proved  difficult.  Several  indivi¬ 
duals  have  made  very  helpful  suggestions  in  this  regard,  however,  and 
further  experimental  work  is  being  carried  out  to  shed  light  on  the 
problem. 

Despite  the  experimental  and  theoretical  difficulties,  the  non¬ 
propagating  hydrodynamic  soliton  is  well  enough  understood  to  be  able 
to  apply  it  co  other  nonlinear  wave  systems.  Future  work  in  this  area 
might  focus  on  the  feasibility  of  an  analogous  acoustic,  light,  or 
flexture  soliton.  Preliminary  work  in  this  area  has  begun  and  may 
be  developed  further. 
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Prediction  of  Tissue  Composition  from 
Ultrasonic  Measurements  and  Mixture  Rules 

Robert  £.  Apfel 

Center  for  Ultrasonics  and  Sonics,  Yale  University, 

?.0.  Box  2159,  New  Haven,  CT  06520 

Abs  tract 

A  methodology  is  presented  for  predicting  the  composition  of  tissues 
from  measurements  of  the  density,  sound  velocity,  and  acoustic  nonlinear 
parameter,  using  mixture  laws  for  the  density,  compressibility,  and  ncniinea 
parameter.  It  is  shown  that  the  mixture  law  for  the  nonlinear  parameter 
plays  an  essential  part  in  this  methodology,  which  leads  to  the  prediction 
of  the  volume  fractions  of  water,  protein,  and  fat  in  a  given  tissue.  Data 
from  the  literature  for  solutions,  blood,  normal  tissue,  and  cancerous 
tissue  are  investigated,  and  predicted  fractions  are  consistent  with  tissue 
compositional  information  available  in  handbooks.  More  experimental  work 
is  needed  with  tissues  of  known  composition  in  order  to  more  fully  test 
the  proposed  methodology. 


RECEIVED: 


1 


Introduction 

-o  what  extent  can  simple  mixture  models  allow  for  the  prediction  of 
the  composition  of  tissues?  For  instance,  simple  mixture  models  for  the 
density,  p,  and  compressibility,  3,  can  be  written  as1 

n 

p  =  2  Prx.  1 

n 

3  =  I  3ixi  2 

i 

where  the  x^  are  the  corresponding  -jclum  fractions  and  n  is  the  number 
of  different  components  that  are  sufficient  ‘o  describe  the  material. 

In  our  tissue  work  we  will  consider  three  basic  components:  i  =  1  is  an  iso¬ 
tonic  water  solution,  i  =  2  is  "pure"  protein  (see  below),  and  i  =  3  is 
pure  fat  (e.g.  lipid  or  oil  base).  To  illustrate  the  utility  of  equations 
1  and  2,  we  consider  the  work  of  Weiser2  and  Weiser  and  Apfel1  who 
measured  density  and  compressibility  of  red  blood  cells  (RBC's)  in  dif¬ 
ferent  saline  solutions,  ranging  from  hypotonic  (74  mM  NaCl,  20  mM  Tris) 
to  hypertonic  (320  mM  NaCl,  20  mM  Tis).  R3C's  can  be  modeled  as  a  mixture 
of  salt  water  and  protein  (i.e.  no  fat).  As  figure  1  illustrates,  the 
linear  mixture  model  works  quite  well.  Moreover,  the  extrapolation  of 
this  data  to  1002  protein  gives  us  values  of  the  effective  density  and 
compressibility  of  pure  protein  as  =  1.37  g/cm5  and  S?  =  0.94x10  11 

cmVdyne.  Note,  also,  that  it  appears  that  the  predicted  sound  velocity, 

1. 

c^,  of  the  mixture,  given  in  terms  of  a*  a. .3  by  c,  =  l/(p.,3„  )', 

does  not  follow  a  linear  mixture  model:  more  data  points,  however,  are 
necessary  to  validate  this  contention. 


ear  parameter, 


Does  this  kind  of  mixture  modeling  work  for  the  nonlin 
3/A,  of  the  solution?  Here  3/A  is  defined  by 

I  ,  ,„c  /dc\  .  d ( 1  / 3  n  _ 

a  !_  0?  Js 

where  (dc/dP)^  gives  the  change  in  sound  velocity  for  an  adiabatic  pres¬ 
sure  change,*  and  d( l/3;/d?  represents  the  material's  ''stiffness"  change 
with  pressure. 

Apfel  recently  derived  an  effective  nonlinear  parameter  of  a  mixture 


a*i(!) .v 


Because  of  S's  dependence  on  the  volume  fractions,  we  see  that  Eq.4  is 
not  linear  in  the  x^'s. 

Equations  1,  2,  and  4  plus  the  condition  that  x^  +  x^  +  =  1 

constitute  the  four  equations  to  be  solved.  Assuming  that  the  tissue 
properties  p,  8,  and  B/A  are  measured,  we  have  four  equations  and  three 
unknowns  —  the  volume  fractions  of  the  aqueous  medium  (water  or  saline 
solution),  protein,  and  fat.  How  does  one  solve  for  the  optimal  choice 
of  the  x^? 


Statement  of  the  Mathematical  Problem 

First  we  repeat  equations  1,  2,  and  4  in  matrix  format,  AX  =  B, 
assuming  the  equations  to  be  exactly  true  and  considering  three  compo 
nents  in  the  mixture, 


Apfel  -  JASA 


* 


S£tB/A)1  3 2 ( 3 / A ) 2  B*(B/A)3 


X3 


3  2  (3/A) 


Later  we  will  find  it  convenient  to  normalize  each  of  the  equations  in 
Eqs.o  by  their  right  hand  sides,  so  that  3  is  a  vector  of  ones.  The  equ 
tiers  in  this  matrix  are  subject  to  the  constraint: 


X1  +  X2  T  X3  =  1‘ 


3ut  the  equations  represented  by  the  matrix  are  only  approximations.  In 
Eq.6,  the  3  vector  components  are  based  on  measured  values  of  p,  3,  and 
3/A.  When  a  set  of  x's  are  chosen,  these  will  produce  effective  pro¬ 
perties,  indicated  by  primes,  that  differ  from  the  measured  values  by 
we  define  the  a's  as  the  normalized  quantities,  as  follows: 


A1  =  ( p  *  — p ) / p 


A2  =  (B'-S)/S 


A3  =  [ (B/A) ’  -  (B/A) ] / (B/A) 


Our  goal  is  to  find  x^,  x2,  x^  such  that  these  errors  are  minimized 


We  choose  the  criterion  that  we  minimize  Q,  where 


Q  =  (A2  +  A2  +  A2). 


We  do  this  by  the  method  of  the  Lagrange  multipier.  The  details,  given 
in  Appendix  A,  result  in  a  set  of  four  equations  and  four  unknowns. 


4 


Input  Parameters 

The  simultaneous  solution  of  the  four  algebraic  equations  or  Eqs.  A-5 
depends  on  the  input  pa: ’meters,  the  elements  in  the  matrix  A,  and  the 
measured  tissue  or  solution  properties,  the  elements  in  the  3  vector. 

Table  1  gives  the  input  parameters  used,  where  T  is  the  temperature  in  °C. 

The  temperature  dependence  of  wafer's  properties  are  taken  from  hand¬ 
books  and,  for  (3/A),  from  the  work  of  Coppens,  et  al . 6  Protein,  con¬ 
sidered  to  be  semi-solid,  is  assumed  to  have  a  negligible  temperature 
dependence  for  density  and  compressibility  over  the  15  Celsius  degree 
span  that  encompasses  all  the  data  to  be  presented.  The  temperature 
dependences  of  the  properties  of  fat  follow  those  of  many  oils.7’*  Namely, 
approximately  -0.001  g/cms/°C  for  the  density  and  about  -3.5  meters/s/°C 
for  the  velocity  of  sound,  from  which  the  compressibility  can  be  calcu¬ 
lated  (a  =  1/pc2).  The  sound  velocity  for  pure  fat  at  20°  is  taken  as 
1451  m/s.  (At  37°C  it  is  estimated  at  1391.5  m/s.) 

Note  that  we  must  assign,  somewhat  arbitrarily,  a  value  for  the  non¬ 
linear  parameter  of  pure  protein.  This  is  true  because  the  weighting 
factor,  3~,  in  Eq.4  is  on  the  order  of  1/20  of  65,  so  that  the  value  of 
3/A  for  the  solution  or  tissue  is  rather  insensitive  to  the  choice  of 
O/A^.  We  choose  a  value  of  10  (at  20°C)  since  many  carbon  based  com¬ 
pounds  have  3/A  values  between  9  and  11. 

For  "pure"  fat,  the  nonlinear  parameter  is  assigned  to  be  11  at  20°C, 
taking  guidance  from  our  earlier  work  in  which  the  nonlinear  parameter 
of  pure  castor  oil  was  measured  at  12.0. 5  The  temperature  dependence  of 
B/A  fat  is  taken  to  be  +.03/°C,  similar  to  that  for  glycerol.5  Our 
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sensitivity  analysis  in  Appendix  3  gives  estimates  of  the  effects  of 
changes  in  our  input  parameters  on  cur  predicted  volume  fractions. 


Results  for  Aaueous  Solutions 


We  new  consider  the  results  based  on  available  incut  data  from  the 


literature. *’ 1 ° ’ 1 1  The  results  for  solutions  are  given  in  Table  2. 


For  dextrose  and  dextran  solutions  we  substitute  for  protein,  su: 
with  ''pure"  density  and  compressibility  given  by  =  1.53  g/cm3  and 


3  *  0.5x10 


cm1 /dyne,  as  determined  from  Eqs. 


1  a r* A  ?  iic  ^ 


using  dens  1  tv  a 


sound  velocity  data  reported  in  the  literature.  O/A)^  is  taken  as  10 
because  the  atomic  make-up  of  proteins  and  sugars  is  similar  in  many 
respects.  (See,  also,  Appendix  3  for  a  sensitivity  analysis.) 

To  convert  the  concentration,  a,  in  grams/cm5  to  volume  fraction, 
we  first  convert  a  to  mass  fraction,  m, 


m  *  a/p,  in  grams/100  grams  solution 


and,  then  using  Eq.l,  the  defining  relation  for  a  mixture  involving 
protein  (or  sugar)  and  water,  we  get 


x^  (protein  or  sugar  volume  fraction) 


x^  (water  volume  fraction) 


10a 


10b 


The  calculated  volume  fractions,  and  those  deduced  from  the  given 
concentrations,  are  expected  to  agree  reasonably  well,  since  p_  and  6„ 


are  deduced  from  the  data.  What  is  perhaps  worthy  of  note,  however,  is 


that  the  calculated  volume  fraction  of  fat  is  negligibly  small,  as 
expected,  and  that  the  measured  value  of  3/A  is,  indeed,  significant  in 
the  calculation.  That  is,  if  B/A  were  7.04  instead  of  6.04  for  dextrose 
at  30°  (a  =  28.9  g/cm3),  the  volume  fraction  of  water  and  protein  would 
each  change  away  from  the  correct  value  by  about  3".  Therefore,  3/A 
results  are  seen  as  part  of  a  consistent  set  of  data  which  determines 


accurate  compositional  information. 

Also  of  note  are  some  negative  volume  fractions,  always  of  small 
magnitude.  These  occur  occasionally  because  the  criterion  we  choose 
for  the  optimal  (Eq.9)  does  not  constrain  the  x,  to  be  positive. 

Results  for  Tissues 

Table  3  presents  results  for  the  components  (x^,  x2>  x^)  of  tissues 
for  which  data  on  density,  sound  velocity,  and  B/A  data  are  available. 
Because  of  the  sensitivity  of  the  results  to  the  accuracy  of  B/A  measure¬ 
ments  (see  Appendix  B),  we  accept  only  B/A  data  having  error  bars  of 
less  than  ±5%.  This  rules  out  all  solid  tissue  B/A  data  employing  the 
finite  element  method. 

There  are  a  number  of  consistencies  in  the  results  that  support  our 
confidence  in  our  methodology.  First,  if  one  looks  at  the  data  for 
multiple  myeloma,’  we  see  changes  of  10%  in  the  B/A  data,  for  a  change  of 
15°C,  whereas  the  corresponding  predicted  change  in  the  volume  fraction 
of  the  predominent  component,  water,  is  only  about  1%.  Since  the  same 
sample  was  used  over  the  entire  temperature  range,  and  since  the  density 
change  is  less  than  1%  over  a  15°  temperature  change,  we  expect  the 
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volume  fractions  to  remain  the  same. 

Second,  the  breast  fat  data  correctly  predicts  minimal  amounts  of 
protein,  and  therefore,  like  the  sugar  solutions  data  which  predicted 
negligible  fat,  conforms  to  reality.  One  anomaly  does  exist  in  the 
measured  data,  namely,  the  sharp  change  in  the  temperature  dependences 
of  the  sound  velocity11  and  nonlinear  parameter.’  Nevertheless,  the  pre¬ 
dicted  tissue  composition  shows  reasonable  consistency  with  temperature. 

Finally,  the  true  test  of  our  approach  is  whether  the  predicted  com¬ 
positions  correspond  quantitatively  to  reality.  Analyses  of  tissue 
composition  are  reported  in  the  literature  and  summarized  in  handbooks. 
Usually  ranges  are  given,  because  the  compositions  depend  on  a  number  of 
factors,  such  as  age,  diet,  state  of  hydration,  etc.  Nevertheless  we  are 
able  to  make  some  comparisons.  For  instance  the  mass  fractions  of  water 
and  protein  for  human  liver  from  the  Handbook  of  Biological  Data11  are 
0.73-0.77  and  0.18,  re^c-er  lively,  with  an  estimate  for  fat  from  animal 
data  being  0.03-0.06.  Our  corresponding  results  (corrected  to  mass 
fractions)  are  0.76,  0.20,  .04.  (The  protein  fraction  is  expected  to  be 
higher,  because  it  includes  other  dissolved  components  that  tend  to 
reduce  the  compressibility.) 

For  cancerous  tissue  it  is  generally  recognized  that  the  water  frac¬ 
tion  increases.  For  instance,  in  rat  hepatoma  the  water  fraction  is 
about  0.71  for  normal  liver  and  0.32  for  the  hepatoma,  whereas  the  fat 
total  is  .062  for  normal  liver  and  0.03 2  for  the  hepatoma  (protein  data  was 
not  given).1*  From  the  data3  we  have  analyzed,  the  water  fraction  goes 
from  0.76  in  normal  human  liver  to  0.90  for  multiple  myeloma,  about  the 
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sane  percent  change  as  for  the  rat  data.  For  fat,  the  mass  fraction  for 
the  human  liver  is  Q.Ou  and  for  myeloma  is  .02.  Thus  the  trends  are  con¬ 
firmed  qualitatively  and,  in  seme  cases,  quantitatively  by  our  predictions. 

General  Discussion 

Should  the  mixture  rules  be  valid?  Or  are  there  substantial  errors 
in  each  cf  them  that  are  in  seme  way  compensated  for  when  the  rules  are 
taken  as  a  set?  v«'e  suggest  that  this  question  depends  on  the  degree  to 
which  molecules  of  the  individual  components  (water,  protein,  or  fat) 
intermix  as  compared  to  "sticxing  to  their  own  kind.”  That  is,  what 
fraction  of  water  molecules  are  surrounded  by  other  water  molecules  and 
what  fraction  of  them  are  in  contact  with  fat  or  protein? 

If  the  intermixing  is  substantial,  as  with  a  2Q7.  mixture  of  methanol 
in  water,  then  the  mixture  rule  is  invalid.  (The  value  of  B/A  for  methanol- 
water  solution  actually  has  a  minimum  which  cannot  be  predicted  by  our 
mixture  model.5)  But  if  large  numbers  of  the  same  atom  groups  are  in 
close  proximity,  as  with  protein  or  fat,  which  have  high  molecular  weights, 
then  the  mixture  rule  is  likely  to  be  valid. 

As  Hartmann15  and  others  have  discussed,  the  nonlinear  parameter  is  tied 
to  the  shape  of  the  intermolecular  or  interatomic  potential  (molecular  scale). 
Thus,  structures  such  as  the  aggregates  of  biological  cells  that  make  up 
tissue,  which  are  much  larger  than  the  intermolecular  scale,  should  not 
play  a  significant  role  in  determining  B/A,  and,  therefore,  B/A  should  not 
provide  information  about  the  hierarchal  structures  of  tissue.  Rather,  by 
providing  an  additional  independent  quantity,  related  to  the  pressure  depen¬ 
dence  of  the  material's  stiffness,  one  is  able  to  characterize  more  fully 
the  composition 
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of  complex  materials,  such  as  tissues. 

If  spatially  dependent  B/A  measurements  can  be  performed  in  vivo, 
as  with  sound  velocity  measurements , 1  *  then  significant  opportunities 
exist  for  diagnosing  pathological  conditions  in  soft  tissues.  Although 
there  have  been  reports  of  such  in  vivo  measurements,15  a  statement, 
made  in  the  description  of  the  work,  namely  that  B/A  is  small  for  fat, 
is  totally  inconsistent  with  vr.  vzivo  measurements,  suggesting  that  the 
nonlinear  parameter  was  not  isolated  from  other  factors. 

Further  work  is  needed  in  accurate  measurements  of  biological  systems 
for  which  the  compositions  are  well  known  and,  therefore,  for  which  com¬ 
parisons  with  our  composition  determination  methodology  are  pcssiole. 
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Appendix  A  Mathematical  Analysis 

We  begin  our  analysis  by  putting  Eq.A  in  a  vector  and  indexed  format 


Q  “  A  •  A 


A-l 


■  ‘Vj  '  Bl>(Aik\  -  V  -  ‘<*1  *  *2  +  *3  -  U 


where  A  is  the  column  matrix  having  A  ,  A^ ,  A^  as  its  components,  and  A  is 
the  transpose  of  a.  Here  the  A's  and  3's  represent  the  elements  in  the  nor¬ 
malized  form  of  Eq.6  (see  Table  1).  I:i  the  format  of  Eq .  A-l  the  double  sub¬ 
script  implies  summation  over  the  repeated  indice,  and  A  is  the  Lagrange 
multiplier. 1 ' 
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'he  minimiza tion  procedure  is  performed  through 


jQ  _  n  alj  -  r  JQ  _  * 

—  -  o,  —  =  0,  - —  -  u 

3X.  aX,  3X-, 

1  4  3 


or 


30 


t—  =  A.  .<$.  (A  ,  x,  -  B  .  )+( A .  .x .  -  B  .  )A  6,  -  X  =  0 

3x  ij  jn  ik  x  i  ij  j  l  ik  kn 


ax . 


where  6 


ja  3x_ 


is  the  Xronecxar  delta. 


Combining  terns  gives: 


30 


■g—  *  A,  (A,,  x,  -  B  )+(A.  X.  -  B .  )A  -  X  =  0 
3x  in  ik  x  i  ij  j  i  in 
n 


2AIn(Aik\  '  V  '  1  *  0 


Including  the  summation  explicitly  gives  the  four  equations: 


2  i*l  Ail(AilXl  +  Ai2X2  +  Ai3X3  ~  V  "  X  =  0 


2  i=l  Ai2(AilXl  +  Ai2X2  +  Ai3X3  "  V  *  0 


2  i=l  Ai3(Ailxl  +  Ai2x2  +  Ai3x3  “  Bi)  A  0 


x  +  x  +  x3  +  (0)X  *  1 


A-' 


A-5 


These  four  equations  can  be  solved  analytically  or  numerically  (see 


Appendix  B)  for  the  four  unknoi/ns ,  x^,  x^,  x^jX,  and  the  A's  C3n  be 


calculated  from  Eqs.  3. 


J. 


Appendix  B  —  Accuracy  and  Sensitivity  Checks 

A  number  of  checks  were  performed  to  test  the  accuracy  and  reason¬ 
ableness  of  cur  formulation: 

1)  Accuracy  of  the  inversion  algorithm 

When  tested  with  simulated  data,  the  algorithm  reproduced  the  cor¬ 
rect  results  within  2  parts  in  10,000. 

2)  Importance  of  3/A  relationship 

If  we  solve  Eq.l  and  2  with  the  constraint  equation,  x.,  +■  x„  -r  x«  = 
for  the  case  of  breast  fat  at  30°C,  we  get  x.  =  -0.096,  x„  =  0.159, 

1  t. 

x  =  0.937,  as  compared  to  the  results  in  Table  3  of  x  =  0.418,  x„  =  .00 
3  1  a 

and  x^  =  .577  when  the  nonlinear  parameter  equation  is  included.  This 
result  points  to  the  necessity  of  using  the  full  complement  of  informa¬ 
tion  available  In  order  to  achieve  reasonable  results. 

3)  Effect  of  input  data 

If  one  or  more  of  the  input  parameters  is  close  to  the  measured 
values  of  p,  8,  or  B/A  for  the  solution  or  tissue,  then  the  results  may 
become  especially  sensitive  to  the  input  parameter.  For  instance,  If 
8^  is  reduced  by  5%,  the  results  for  breast  fat  at  22°C  change  as  follows 
water  fraction  reduced  by  19%,  fat  fraction  up  16%,  and  a  significant 
change  in  the  protein  fraction.  Similar  changes  occur  for  measured  B/A 
values  that  approach  the  value  of  B/A  for  pure  fat  (-11).  The  effect 
is  smaller  for  tissues  for  which  the  B/A  value  is  lower,  as  with  most 
other  tissues  (e.g.  liver).  But  the  effect  is  not  negligible.  For 
instance,  if  we  take  the  liver  data  and  reduce  the  measured  density  by 
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10r: ,  the  predicted  water,  protein,  and  fat  volume  fractions  change  by  0.3%, 
2.9%,  and  3.5%  respectively.  The  corresponding  changes  for  a  10%  increase 
in  3/A  are  7.5%,  10.6%,  and  over  100%. 

This  sensitivity  check  reinforces  our  caveat  that  accurate  m  sure- 
ments  of  tissue  properties,  such  as  are  possible  with  phase  methods  for 
3/A  ,  are  a  prerequisite  for  getting  reliable  compositional  infor¬ 

mation.  Moreover,  we  should  measure  3/A  of  lipid  oils  in  order  to  get 
a  more  accurate  estimate  of  the  nonlinear  parameter  of  pure  fat,  thereby 
reducing  cur  uncertainties  in  tissue  composition  especially  for  those 
tissues  containing  a  significant  fat  component. 
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TABLE  1  -  INPUT  PARAMETERS 


3  A, ,  =  P . 
11  I 


•  AU  '  °2  ’ 


I  1.005  -  .00029*(T-2Q)  for  tissues  and  whole  blood 
|  .593  -  .0CG29x(T-2Q)  for  dextrose,  dextran,  and  3SA  sciatic: 


1.37  for  "pure"  protein 

1.53  for  "pure"  sugar  (for  dextrose  and  dextran  solutions) 


3  A,  =  P,  =  0.33  -  0.001*(T-20) 

x  J  j 

1.513*10s  t  24G*(T-30);  saline 

-  A„ ,  *  3,  =  1 ( p  c.*),  where  c1  = 

“*  x  11  1  1.510x10s  -  240x(.T-30);  sugar 

0.94  for  pure  protein 

3  \  =s  3  = 

2 

0.5  for  pure  sugar  (for  dextrose  and  dextran  solutions) 
a  a^3  =  03  =  5.4(1  +  .00463x(T-20)>,  where  c3  =  1.451x10s  at  20°C 

S1(3/A)A31  =  312(3/A)1,  where  (B/A^  *  5.3  +  0.02x(t-30) 

Bz(B/A)A32  =  B^2 (B/A)^ ,  where  (B/A)2  =  10.0  +  .031x(t-20)  for  protein  and  sugar 
a2(3/A)A„  =  B,2(B/A)3,  where  (B/A),  =  11.0  +  0.031x(T-20) 


Densities  in  g/cm3;  compressibilities  in  cm2/dynexl0 1 1 ;  T  in  C;  c  in  cm/s. 
Note  that  all  A.'s  are  normalized  so  that  all  B,'s  are  1. 

r 

This  value  is  consistent  with  compressibilities  of  solid  polymers,  such 
nylon,  polystyrene  and  lucite. 


This  value  follows  from  the  compressibility  of  solid  dextrose,  as  measured 
by  Bridgman. 1 5 


